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For many real spin-glass materials, the Edwards- Anderson model with continuous-symmetry spins 
is more realistic than the rather better understood Ising variant. In principle, the nature of an oc- 
curring spin-glass phase in such systems might be inferred from an analysis of the zero-temperature 
properties. Unfortunately, with few exceptions, the problem of finding ground-state configurations is 
a non-polynomial problem computationally, such that efficient approximation algorithms are called 
for. Here, we employ the recently developed genetic embedded matching (GEM) heuristic to investi- 
gate the nature of the zero-temperature phase of the bimodal XY spin glass in two dimensions. We 
analyze bulk properties such as the asymptotic ground-state energy and the phase diagram of dis- 
order strength vs. disorder concentration. For the case of a symmetric distribution of ferromagnetic 
and antiferromagnetic bonds, we find that the ground state of the model is unique up to a global 
0(2) rotation of the spins. In particular, there are no extensive degeneracies in this model. The 
main focus of this work is on an investigation of the excitation spectrum as probed by changing the 
boundary conditions. Using appropriate finite-size scaling techniques, we consistently determine the 
stifi^ness of spin and chiral domain walls and the corresponding fractal dimensions. Most noteworthy, 
we find that the spin and chiral channels are characterized by two distinct stiffness exponents and, 
consequently, the system displays spin-chirality decoupling at large length scales. Results for the 
overlap distribution do not support the possibility of a multitude of thermodynamic pure states. 

PACS numbers: 75.50.Lk, 64.60.Fr, 02.60.Pn 



I. INTRODUCTION 



The discovery of materials without magnetic ordering 
down to zero temperature, but with an unusual cusp in 
the non-linear magnetic susceptibilit}^ indicating a tran- 
sition to a magnetically frozen but random state, has 
prompted a theoretical effort spanning over thirty years 
and aimed at the understanding of the perplexing prop- 
erties of these spin-glass systems^i^i^i^. Albeit mostly 
driven by the theoreticians' hope to understand the pe- 
culiar behavior of such systems within the framework of 
statistical mechanics, the results of this research have 
found widespread application in seemingly rather distant 
fields such as the theory of associative memory^, mod- 
els of the immune system! or error correcting codes-'-. 
Spin-glass science thus provides for a classical example 
of how research driven by curiosity rather than applica- 
tion can lead to unexpected progress in diverse research 
areas. The canonical case of a spin glass, i.e., a noble 
metal host doped with magnetic transition metal ions^ 
is not of particular interest as a material apart from the 
peculiarity of its properties. In recent years, however, 
many new materials with possible spin-glass phases have 
been discovered in compounds with highly frustrating 
lattice structures''^. The large degeneracy of the ground 
state arising from geometric frustration alone gives rise 
to a host of novel and exotic magnetic phases, including 
spin glassesii, but also spin liquids'^'^^ and spin icesi^. 



Experimental evidence suggests that spin-glass like phe- 
nomena appear to lie at the heart of several (potentially) 
technologically important phenomena, such as high tem- 
perature superconductivity, colossal magnetoresistance, 
and the anomalous Hall effect. 

Other than might be suspected from the pure vol- 
ume of published research in the field-'^, however, 
owed to the unwieldy nature of the problem, even 
some very basic questions such as the existence of a 
finite-temperature phase transition have in general re- 
mained unanswered to date^. Starting from the, by 
now, rather complete'- understanding of the unusu- 
ally complex mean-field theor y'^1'^1'^ of the fully con- 
nected (infinite dimensional) Sherrington-Kirkpatricki^ 
(SK) model, recent research has focused on the case 
of finite-dimensional systems. For the nearest-neighbor 
0(n) Edwards- Anderson (EA) modeP° with Hamiltonian 

H = -^Jij Si ■ Sj, (1) 

which has been studied in most detail, most of the effort 
has been devoted to the discrete case of classical Ising 
spins (n = 1). Even for this simplest case some controver- 
sial questions, such as the existence of a glassy phase at 
finite-temperatures in three dimensions (3D), could only 
be definitely answered in very recent years'^. In general, 
many features such as the generic structure of valleys and 
barriers in the free energy and, more generally, the ques- 
tion of whether spin glasses in finite dimensions behave 
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essentially alike or differently from the mean-field case 
remain poorly understood^ii^. Even the allegedly sim- 
plest case of the Ising spin glass in two dimensions (2D) 
remains a topic of lively research activity2^i2ii^i2£i22. 

For most of the materials featuring glassy magnetic 
phases, however, continuous XY or Heisenberg spins cer- 
tainly are a more accurate description than the Ising 
model of extreme anisotropy in spin space^i^i. From 
the perspective of condensed matter physics, much at- 
tention has been paid to the experimentally important 
case of weak, but frustrating impurities and the question 
of a resulting loss of coUinear order. Isolated frustrat- 
ing antiferromagnetic defects in an otherwise ferromag- 
netic system couple to each other via an effective dipolar 
interaction22i2^j but it remains an open question whether 
collinear order can be destroyed with an arbitrarily small 
concentration of sufficiently strong defect bonds, and 
what role is played by the spatial dimension'^" i^i^^i^. 
Such models are potentially relevant to the weakly doped 
copper-oxide parents of high-Tc superconductors^^i^iiS^, 
while periodically frustrated XY systems are models of 
Josephson junction arrays^ii^i^li^. For the strong dis- 
order case of a spin glass with symmetric bond distribu- 
tion, the question of whether there is a non-zero tran- 
sition temperature (for either or both of the spin and 
chiral sectors) in 3D XY and Heisenberg EA models is 
currently (again) hotly debate d'^^'^^i'^^'^^i'^^'^'^ . 

For the SK model, the ordered phase is characterized 
by a breaking of the abstract replica symmetry^, but the 
nature of symmetry breaking in the short range case is 
less clear. Starting out from an essential analogy with the 
behavior of the ferromagnet, a phenomenological droplet 
scaling theory for the effective coupling in spin glasses 
has been formulated^i^'^i^'^. The surface (free) en- 
ergy of droplet excitations of length L is assumed to scale 
as AE ~ defining the spin stiffness exponent 6s- 
Put into the context of a renormalization group (RG) 
frameworl^, 9s determines whether a system scales to 
weak coupling {Og < 0) such that droplets can always 
destabilize the ordered phase, thus restricting its exten- 
sion to zero temperature, or scales towards strong cou- 
pling {9s > 0) such that large-scale excitations are suf- 
ficiently costly to allow for a stable spin-glass phase at 
finite temperatures < T < Tsg- In the infinite-range 
SK model, the notion of droplets is maldefined and, in- 
stead, it is found that excitations of arbitrary size can 
be invoked at constant energyi^iiiiiSii^. Consequently, if 
mean-field theory applies to low-dimensional spin glasses, 
low-energy excitations are of rather different nature im- 
plying, for instance, space filling domain walls^^. Hence, 
an analysis of overlaps between ground and systemati- 
cally excited states can in principle allow to distinguish 
between the droplet and mean-field picture o^^'^^ . A stan- 
dard way of producing excitations is a change of bound- 
ary conditions (BCs), e.g. from periodic to antiperiodic, 
the energy difference between both configurations cap- 
turing the energy of a relative domain wall^i^. 

For the 2D Ising spin glass, domain- wall energies scale 



with 6s = —0.287(4) for Gaussian bond distribution— 
and ds — for bimodal couplings^, and a true thermo- 
dynamic spin-glass phase occurs only at zero temperature 
in both cases. From the RG ansatz of Ref. deter- 
mines the divergence of the correlation length ^ as T ^ 
via v — —l/9s, which is found to be consistent with 
finite-temperature Monte Carlo results for the Gaussian 
case^^i^. For the bimodal ± J distribution, on the other 
hand, the formal v = oo has been interpreted as exponen- 
tial divergences^ of ^, but some evidence for an algebraic 
singularity has also been presented recently2£ (see also 
Ref. IstI ). For the case of continuous XY or Heisenberg 
spins, the situation is complicated by the fact that the 
non-collinear spin structure in the ordered phase allows 
for the distinction of proper and improper 0(n) transfor- 
mations, such that the continuous SO(n) rotational sym- 
metry is augmented by a Z2 or Ising-like chiral symmetry 
related to the determinant of the 0(n) rotation matrix^-. 
It has been argued!^ that these chiral degrees-of-freedom 
might decouple from the rotational ordering, resulting in 
different transition temperatures (above the lower criti- 
cal dimension) or, at least, an additional chiral stiffness 
exponent 9c (for Tsg = 0)- Such a scenario is — after 
a long-winded debate — now quite well established in 
some periodically frustrated magnet o^°i^^'^^ , but remains 
fervently debated for the case of spin glasses^SiildSiiiiiS. 
For the 2D XY spin glass, neither has dg been determined 
consistently, with estimates ranging from^^ 9s = —1.0 
tc^^ 9s = —0.4, nor has the question of a possible spin- 
chirality decoupling been satisfactorily resolved^^'^^. 



Methodologically, the above uncertainties regarding 
the XY spin glass are related to various technical diffi- 
culties. Firstly, finding ground states of spin-glass model 
systems is a computationally hard problem, generically 
believed to require an effort growing exponentially with 
system size^. The 2D Ising spin glass on planar graphs 
is the only non-trivial exception to this rule, being poly- 
nomial computationally^*, and allowing for the rather 
precise estimate of 9s cited above. Thus, reliable spin- 
glass ground state computations depend in general on 
rather elaborate and specifically tailored approximation 
algorithms to treat reasonably sized systems. In this 
context, we use here the recently proposed "genetic em- 
bedded matching" (GEM) technique^i^Si to investigate 
the zero-temperature properties of the 2D XY spin-glass 
model. For the Ising problem mentioned, due to the exis- 
tence of rather pronounced finite-size effects^i relatively 
large system sizes turned out to be necessary for reliable 
estimates, e.g., of 9s. Below, we will see how progress 
in accessing much larger system sizes using the GEM 
approach now allows for investigations where finite-size 
effects seem for the first time rather well controlled. Sec- 
ondly, there is no full agreement as to how BCs should 
be chosen in order to precisely excite either spin or chiral 
defects^!^^. This is of particular concern here, given that 
strong corrections to scaling depending on the chosen set 
of BCs have been observed for the Ising cas o^'*'^'^ . Here, 
we systematically study a wide range of BCs and resolve 
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the boundary dependent corrections to scaling. 

The rest of the paper is organized as follows. Section 
In] briefly discusses the employed GEM technique and the 
technical set-up of the calculations. In Sec. lIIIl we discuss 
the asymptotic ground-state energy and the spectrum of 
metastable states found in the GEM calculations, while 
Sec. IIVI reports on the properties of excitations induced 
by a change of BCs. In particular, spin and chiral stiff- 
ness exponents and the associated fractal dimensions of 
domain walls are determined. Some of these results al- 
ready were reported in part in Rcf. 6^. Section |V] con- 
tains our results for the scaling of overlap distributions, 
and Sec. |Vl] presents some results on the distribution of 
internal fields and the disorder strength vs. disorder con- 
centration phase diagram of the 2D randomly frustrated 
XY magnet. Finally, Sec. IVIII contains our conclusions. 

II. ALGORITHM AND TECHNICAL SET-UP 

Zero-temperature properties might be inferred from 
the T — > limit of finite-temperature calculations such 
as Monte Carlo simulations. The latter, however, suffer 
from severe slowing down as the T = critical point is 
approached for this system, thus rendering the extrapo- 
lated results rather unreliable. Hence we chose, instead, 
to determine ground states directly from non-equilibrium 
methods. Apart from the 2D Ising EA model on pla- 
nar graphs, such computations are believed (and in some 
cases proven) to be NP hard^^^ optimization problems^^. 
To make progress with the 2D XY model, we take inspi- 
ration from the approach making the 2D Ising EA model 
solvable^^. In particular, we embed Ising variables into 
the continuous XY spins: with respect to a random, but 
fixed direction r in spin space, the 0(n) spins Si of ([J) 
decompose as Si = s} -I- S^ = {Si ■ r)r + S^. This 
results in a decomposition of 7i as 7i = Ti'"'" +TL^'^ with 

^'■'"=-E4<4' (2) 

where e[ — sign(S'i • r), and where the effective couplings 
Jlj are given by 

4 = J.,|5,.r||5,.r|. (3) 

If now the original 0(n) rotations of the spins are re- 
stricted to reflections along the plane defined by r, J^j as 
well as are easily seen to be invariant under such re- 
fiections of individual spins, whereas the change sign. 
Hence, for a given r, Ti'"'" in ([2]) takes on the form of an 
Ising model with effective Ising "spins" 

The tractability of the 2D Ising ground-state problem 
in polynomial time comes about via the transformation 
to a graph-theoretical formulation=S* . It rests on a dual- 
ization of the model: define elementary plaquettes of the 
square lattices to be frustrated if one or three of their 
surrounding bonds are antiferromagnetic-^. For a given 



spin configuration, mark dual bonds of the square lattice 
as "active" if the corresponding original bond is unsatis- 
fied, i.e., if JijSiSj < 0. Then, from the very definition 
of a frustrated plaquette, the set of "active" (dual) bonds 
consists of individual chains which are either closed loops 
or open chains connecting pairs of frustrated plaquettes, 
and the total energy oi TC in ^ is (up to a constant) 
identical to the total sum of along the "active" bonds. 
Since on a planar lattice closed loops can always be re- 
moved by contraction, a ground state of the spin sys- 
tem corresponds to a minimum- weight perfect matching 
of the frustrated plaquettes, i.e., a configuration of the ac- 
tive bonds matching up all pairs of frustrated plaquettes 
while minimizing the total weight ^ of bonds con- 
tained. Such problems can be solved in polynomial time 
by means of Edmonds' "blossom algorithm" and, most 
importantly, the solution can be shown to always corre- 
spond to a valid spin configuration for the case of planar 
graphs^'*. As a result, application of (anti-)periodic BCs 
is restricted to at most one direction to preserve pla- 
narity. See Refs. [stIIssIItoI for a more detailed account of 
the relation between the Ising ground state and graph- 
theoretical matching problems. 

Using the decomposition ^ for the 0(n) Hamiltonian 
([1]), the blossom algorithm can be used to find exact 
ground states of the embedded Ising spins e^, correspond- 
ing to a reflection of some of the continuous spins with re- 
spect to the plane defined by a randomly chosen direction 
r. A random sequential sampling of directions r is used 
to statistically recover the full 0(n) symmetry. This em- 
bedded matching procedure is a strictly downhill search 
which, however, due to its non-local nature originating 
from the exact ground-state computation for the embed- 
ded Ising spins, features far less metastable states than 
a zero-temperature quench with any local dynamics^^. 
Due to the dependence of the effective couplings J^j of 
Q on the configuration {Si}, the embedded matching 
dynamics normally drives the system into a (low-lying) 
metastable state instead of the true ground state. 

To find true ground states, the above embedded match- 
ing must be supplemented by a complementary protocol. 
We do so by inserting the embedded matching (i.e., mini- 
mization technique) as an optimization component (sub- 
routine) into a specially tailored genetic algorithm§^2^. 
Here, a whole population of candidate ground state con- 
figurations is being worked on, continuously undergo- 
ing an improvement process mimicking natural evolution, 
which allows for a very efficient exploration of the land- 
scape of metastable states of the system: 

• Crossover: a random pair of spin configurations 
is selected as "parent configurations" and locally 
rigid domains or clusters are exchanged at ran- 
dom between them to produce "offspring" (see also 
Sec. HUB]). 

• Mutation: a small proportion of spins of either 
offspring are randomly perturbed. 
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• Optimization: both offspring are optimized using 
the embedded matching technique. 

• Selection: each parent is replaced by the morpho- 
logicaUy closer offspring in case the latter has a 
lower internal energy. 

The most crucial step in this sequence is the "crossover" 
operation. Indeed, special attention must be paid to 
the way configurations are combined in the crossover to 
balance the need for "genetic" diversification with the 
preservation of the high optimization on short length 
scales already achieved at intermediate steps of the pro- 
cess. Statistical analyzes show that the above combi- 
nation of techniques in this genetic embedded matching 
(GEM) approach allows for a reliable determination of 
ground states for 2D XY spin glass samples of up to 
about 30 X 30 spins with current computing resourcesi^^. 
The algorithmic details of the GEM technique are quite 
intricate. A presentation of those details here would de- 
tract from our main agenda, which is to focus, and report, 
on the physical properties of the ground state of the 2D 
XY spin glass. A full account of the algorithmic details 
of the GEM technique can be found in Ref. [tO- 

For specificity, we assume a bimodal distribution of the 
quenched random couplings Jij of ((T|), 

P( J,,) ^{l-x) 6{J,j -J) + x 5{J,j + AJ), (4) 

where A denotes the relative strength of the antiferromag- 
netic impurity bonds Jij = — AJ and x their concentra- 
tion. Unless stated otherwise, we averaged with respect 
to the distribution (U by taking 5 000 independent disor- 
der realizations and concentrated on the symmetric point 
X — 0.5, A = 1, deep in the spin-glass regime of model 
([1]) with P{Jij) from Eq. ([4]). All computations discussed 
here were performed on square lattices. Choosing the pa- 
rameters of the GEM runs to arrive at true ground states 
with high reliability, the run-times for single ground-state 
computations on a single 2.8 GHz PentiumlV CPU are 
between a few seconds for the smallest lattice sizes up to 
about 8 hours for a 28 x 28 system. In total, production 
of the data presented here consumed the equivalent of 
about 80 years of single CPU time. 

III. ENERGETIC PROPERTIES 

A. Asymptotic ground-state energy 

To determine the asymptotic ground-state energy, we 
performed a finite-size scaling analysis with systems of 
different BCs, namely open and open-periodic (as men- 
tioned above, fully periodic boundaries are not accessible 
to the method due to the planarity constraint^). We 
used L X L systems with linear dimensions L = 6, 8, 10, 
12, 16, 20, 24, and 28 and with 5 000 disorder replica 
for each system size. We consider the average, size- 
dependent internal energy per spin, e(L) = {E/Nq)j, 



where denotes a disorder average and Nq is the 
number of lattice sites. From generic finite-size scaling 
arguments2^, at a T = critical point, e(L) should scale 

79 

e{L) - Coo - L-(^-''\ (5) 

where d denotes the spatial dimension and 9s = —l/v 
was assumed^. As usual, corrections to this form can be 
either non-analytic, stemming from irrelevant operators, 
or analytic and proportional toL^*^, fc=l, 2, from 
non-linear scaling fields*". From the e-expansion of the 
Ising spin glass, the leading non-analytic correction is 
(to first order in e = 6 — d) found^^ to be proportional 
to L^'~^, asymptotically clearly smaller than the leading 
analytic correction term. Assuming a similar behavior for 
continuous spins, we shall restrict ourselves to analytical 
corrections here. Then, the most general form would be 

e(L)-eoo = AL^'^''-'^^'>+BL-^+CL-^+DL-^+. ... (6) 

We first consider the case of fully open boundaries. Not- 
ing that in this case the analytical corrections can be at 
least partly attributed to the presence of edges and cor- 
ners (cf., e.g., Ref. Isil) . Campbell et alr^ take a slightly 
different approach and, instead of using e(L), consider 
the energy e^'^{L)/2 per bond as the fundamental prop- 
erty, making the following scaling ansatz, 

e'^''{L)-e^J = B*/L + C*/L', (7) 

with B* representing edge effects and C* accounting for 
corner corrections. No term proportional to L~'^'^~^=) 
is included in Eq. ([7]) arguing that free boundaries do 
not force any domain walls into the system. Indeed, ex- 
cellent agreement with the data is found under this as- 
sumption for the case of the Ising spin glass^'-. Since 
Le{L) = {L - 1) e'^'^(i), the parameters of Eqs. (O and 
O are related as Coo = e^, B = - ej^, C = C* - B* , 
D = —C*, and A — 0. Note that this effectively reduces 
the number of fit parameters of the form ([5]) with A ~ 
from four to three, the 1/L'^ term merely being produced 
from the 1/L^ correction of ([7]) by the transformation 
from bonds to sites. A fit of the form ([6|) with A — 
to our data for free boundaries yields Coo — —1.5520(14) 
with good quality-of-fit Q = 0.35. Using B* and C* as 
free parameters, on the other hand, yields a very similar 
estimate Coo = -1.55 331(56) with Q = 0.48. This last 
fit is shown together with the data in the upper panel of 
Fig. [1] It is interesting to note that the resulting param- 
eter estimates B* = -0.375(17) and C* = 0.059(98) are 
in fact rather close to the values B* ~ —0.3205(9) and 
C* — 0.042(3) found for the Gaussian 2D Ising case and 
much further away from the estimates B* = —0.5492(20) 
and C* = 0.506(18) for the bimodal Ising spin glasses. 
This is consistent with the expectation that the continu- 
ity of the XY spins effectively smoothes out the discrete 
bimodal distribution of couplings Jij as far as energetic 
defects are concerned. 
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FIG. 1: (Color online) Top: Average ground-state energy per 
site for open-open and open-periodic boundaries. The fits 
are described in the main text, and the horizontal line shows 
the estimated asymptotic value of Coo = 1.55 331. Bottom: 
domain- wall related correction to the internal energy for open- 
periodic BCs, which we denote by eDw{L), and which is given 
by eow{L) = e(L) - e^^ - {B* - e°^)l{2L) + B*/{4L^). Here, 
the open-boundary estimates Boo = 1-55 331 and B* = —0.375 
have been used. The line shows a fit of the form euw{L) — 
Al-(d-e,)^ yielding an estimate 9^ = -0.432(94). 

For (mixed) open-periodic boundaries there is no rea- 
son to exclude an term related to domain-wall 
trapping. In the scaling ansatz ([7|) for the bond ener- 
gies, the correction should be omitted since the 
system has no corners, and the edge contribution should 
be cut in half. Since now 2Le{L) ^ {2L - l)e'=^(i), 
one identifies B = {B* - e^)/2 and C = -B*/4:. The 
statistical precision of the data is not sufficient to per- 
form an unrestricted fit of the general form ([6]) . Fixing 
C = D — with the four unrestricted parameters Coo, 
A, 0g and B remaining, we arrive at e^a — —1.5525(13), 
9s = -0.49(69) and Q = 0.35. Fixing only D = 0, but 
reducing the number of parameters to four again by us- 
ing the parametrization with B* and C* yields a rather 
similar estimate Coo = —1.5529(13), 9s = —0.59(77) and 
Q = 0.40. Using the estimates Coo = —1.55 331(56) and 
B* = —0.375(17) from the open-boundary case to reduce 
the number of parameters to two, on the other hand, 
yields a good fit {Q ~ 0.38) with a much more precise es- 
timate 9s — —0.432(94) (note, however, that the quoted 
statistical error does not take the uncertainty in the fixed 
parameters Coo and B* into account) . This fit is displayed 
together with the data in the upper panel of Fig. [TJ The 
lower panel shows the domain- wall related part eDw(-^) 



of the energy correction for the open-periodic case. As 
is clearly seen, the domain-wall contribution is a very 
small correction to the internal energy. The estimate for 
9s is compatible with the result found from the domain- 
wall calculations for periodic BCs at fixed aspect-ratio 
R — 1 presented below in Sec. IIV[ but too large in mod- 
ulus compared to our final estimate from the aspect-ratio 
scaling. One interesting question is whether spin and/or 
chiral defects should contribute to the domain-wall cor- 
rection to the internal energy. We will show below in 
Sec. IIVI that the chiral defects are of much lower abso- 
lute energy for the considered lattice sizes than the spin 
defects, such that essentially only the spin defects are at 
play here, justifying post factum the interpretation of 9s 
as the spin-stiffness exponent. Our best estimate for the 
asymptotic ground-state energy, eoo = —1.55 331(56), 
should be compared with the value Coo = —1.401 97(2) of 
the bimodal 2D Ising spin giass^ which is, perhaps sur- 
prisingly, only about 10% higher. Hence, the enhanced 
ability of the continuous model compared to the discrete 
Ising case to locally adapt does not allow to substantially 
relieve the frustration imposed by the bond disorder (the 
unfrustrated system yielding, of course, eoo = —2). 



B. Spectrum of metastable states 

While a finely-tuned choice of parameters in the GEM 
approach guarantees high reliability for the optimization 
to arrive at a ground state, less optimal settings (in par- 
ticular a reduction of the initial population size below 
a certain threshold) result in runs converging to low- 
lying metastable states instea d^^'^" . Additionally, and 
irrespective of the parameter choice, a part of the spec- 
trum of low-lying metastable states appears as interme- 
diate population members in the course of the optimiza- 
tion procedure. For a finite system this spectrum ap- 
pears discreteiii (see below) , and one naturally wonders 
what the nature of these excitations is. Comparing the 
configurations of some of these states, it is clearly seen 
that on short length scales the system is rigid within 
the space of energy minima, i.e., it naturally decomposes 
into stiff clusters of solidary spins whose individual mo- 
ments have very similar orientations in all of the low- 
lying metastable states, apart from relative, essentially 
exact 0(n) rotations of whole clusters. Indications of 
such rigid clusters have been found in previous studies 
of IsingSiS and continuous-spin system o^^i^^'^'* . How- 
ever, the genuine existence of such clusters for continuous 
spins had so far not been as well established numerically 
as we are able to achieve here with the GEM method. 
Note that such clusters are the finite-energy analogue 
of the "backbone" or "rigid lattice" of sites with con- 
stant relative spin orientation in all ground states ob- 
served in the ±J Ising spin glass with a large ground- 
state degeneracj*^. The GEM approach, and ultimately, 
its ability to determine the ground states, heavily relies 
on this cluster structure in the choice of operation for 
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FIG. 2: (Color online) Local rotation matrices between the 
ground and first excited states for two 28 x 28 disorder realiza- 
tions with open-periodic boundaries (periodic boundaries in 
the vertical direction). The arrows indicate the rotation an- 
gles, and blue squares are drawn in the case of improper rota- 
tions, i.e., negative determinant of the 0(2) matrix. See below 
in Sec. |V] for details about how these matrices are computed. 
The relative excitation energies are I^EjE « 1.5 X 10"^ (up- 



per configuration) and A.E/E ; 
tion), respectively. 



2.4 X 10 (lower configura- 



the crossover procedure exchanging such rigid domains 
between the "parent" replica, cf. the discussion in Sec. |TT] 
and Ref.^'3,. As a result, the GEM algorithm effectively 
operates directly on the space of metastable states in- 
stead of on the v^fhole of the continuous phase spaceS^. 
Products of independent 0(n) transformations of adja- 
cent clusters (plus local excitations close to the bound- 
aries between clusters) correspond to the domain-wall 
or droplet excitations postulated by the droplet theory 
of the spin glass phase^i^i^ and directly observed in 
Ising spin glasses^Si^iZl. Finally, from the perspective 
of defect theory, the continuous symmetry of the order 
parameter also allows for topological excitations, here, 
according to the homotopy groups of 0(2), reahzed as 
domain-walls and vortex excitations^^. These can be ex- 
plicitly observed in configurations^'''^'^, contributing to 
the thus rather rich variety of excitations appearing in 
this system. 

Although a multitude of states with energies just 



slightly above the ground state (with, e.g., AE/E as 
low as « 10~^ for some 16 x 16 systems) is found, these 
energy differences are still far above the numerical reso- 
lution limits, and the configurations can be shown for 
each case to be truly different (see, for instance, the 
state differences depicted in Fig. [2|) . Naturally, each of 
these states has an |0(n)| degeneracy corresponding to 
a global (proper or improper) rotation of all spins. Con- 
versely, however, we checked for additional non-trivial de- 
generacies and found that statistically independent GEM 
runs always converged to an 0(n) rotated copy of the 
same state. Hence, the ground state of the system ap- 
pears unique modulo the global symmetry. Consequently, 
quantities related to statistical fluctuations vanish at 
T = and rj, the exponent of the correlation function, 
is equal to zerc^^. In terms of the clusters of rigid spins 
observed, the unique ground state results from the fact 
that the rigid clusters arc not free, but their rigid 0{n) ro- 
tations have some non-vanishing excitation energy from 
necessary local adaptations of boundary spins. 

By design, the GEM approach does not explicitly find 
all states of a given energy. In particular, as a non- 
equilibrium technique it does not find states with prob- 
abilities corresponding to their respective Boltzmann 
weights. However, there is no reason to believe that 
in all cases a single of possibly a multitude of degen- 
erate ground states would dominate so strongly that not 
a single degenerate state appears in the thousands of in- 
dependent runs undertaken. In contrast, the combined 
binary symmetry of the bimodal Ising spin glass results 
in a massive degeneracy of the ground stated, such that 
already a typical 10 x 10 system has between 10"* and 
10^ ground states®^. Monte Carlo simulationsS^iSS and 
Migdal-Kadanoff RG calculationsSS have led to specu- 
lations that the pairing of bimodal couplings with the 
discrete chiral symmetry might give rise to a similar sit- 
uation in the bimodal XY spin-glass model considered 
here. Certainly, if calculations such as finite-temperature 
Monte Carlo simulations are sensitive to a (possibly nar- 
row) range of energies, the observed dense spectrum of 
states might effectively appear as quasi-degenerate, re- 
sulting in an erroneously non-zero estimate of 77 from 
such simulations^iS^. Also, one might speculate about 
degeneracies in excited states which we have not inves- 
tigated here, but might have shown up in some of the 
previous studies. 

The low-lying metastable states appearing in the pro- 
cess of a GEM run might allow to learn something about 
the spectral properties of the system. Assuming a unique 
leading exponent in the scaling of defect energies with 
length scale, it is interesting to see how the lowest- 
energy excitations corresponding to the gap relate to the 
domain-wall or droplet excitations that are crucial for 
the understanding of the spin-glass phase. In Fig. [3] 
we show the gap energies for a series of L x 6L sys- 
tems with open-periodic BCs (the non-trivial aspect ra- 
tio is chosen to reduce corrections to scaling, see the 
discussion below in Sec. lIVp . A fit of the functional 
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form (Ai?gap)j ^ L^ffp to the data for L > 5 yields 
an estimate of 6'gap — —1.447(31). Thus, the energy 
gap closes rather rapidly, leading to the expected quasi- 
continuous spectrum of metastable states in the thermo- 
dynamic limit. Also, 0gap is clearly far more negative 
than the estimate 6s ~ —0.3 of the spin-stiffness expo- 
nent found below in Sec. IIVI from the domain- wall energy 
scaling. In fact, the energies of domain walls induced by 
a change from periodic to antiperiodic boundaries, addi- 
tionally shown in Fig. [3] for comparison, are much larger 
than the gap energies for the whole range of system sizes 
considered. It is thus clear (at least as far as the range 
of available finite lattice sizes allows to judge) that the 
excitations reflected in the gap scaling are not in general 
extended defects of the domain-wall type. If, as argued 
here, metastable states are related to the ground state by 
rigid 0{n) rotations of one or several clusters (plus cer- 
tain local adaptations at the cluster boundaries), it ap- 
pears that the gapped excitations correspond to rotations 
of single (or few) such clusters, resulting in a lower excess 
energy than observed from the larger number of clusters 
that recombine to form a domain- wall defect. This view 
is corroborated by the snapshots of defects corresponding 
to the differences between ground and first excited states 
depicted in Fig. (5] Besides the extended defects with sizes 
of order L which might as well be incurred by a change 
of BCs (upper panel in Fig. [5]), in many cases much more 
localized defects, corresponding to single rigidly 0(n) ro- 
tated clusters, are the lowest excitations (lower panel in 
Fig. (2). These two situations (represented by the upper 
and lower panels of Fig [2]) allow us to interpret the rapid 
drop of gap energies indicated by the estimate of 0gap as 
the effect of an additional implicit condition of minimum 
energy. For a given disorder configuration, Ai?gap is not 
the average energy cost incurred by the 0(n) rotation 
of a single cluster, but the minimum cost of all possible 
such rotations of clusters. Finally, it should be noted 
that AE'gap is naturally not confined to either spin or 
chiral^i^ sectors, such that depending on the disorder 
configuration spin or chiral defects make up the gapped 
excitation. 



IV. DOMAIN WALLS 



For a systematic determination of the spin and chiral 
stiffness exponents 9s and 9c, respectively, a controlled 
and direct insertion of the corresponding defects into the 
system is needed. This is most conveniently achieved by 
using an appropriate choice of BCs. Droplet excitations 
could be created by similar means^liSS^ but we have not 
explored this yet. 
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FIG. 3: (Color online) Scaling of the energy gap between 
ground state and lowest excited metastable state as esti- 
mated from the GEM approach for Lx6L systems with open- 
periodic BCs. The line shows a fit of the form {Ai?gap)j ~ 
j^^gap to the data. For comparison, the upper data set shows 
the average energy of domain walls induced by a change from 
periodic to antiperiodic BCs. 



A. Choice and implementation of boundary 
conditions 

For the Ising spin glass, defect energies are carried by 
domain walls or droplet surfaces which are the sharply 
localized boundaries of regions of inverted spins. For 
continuous symmetry, the notion of defects is more sub- 
tle since the spin direction varies smoothly, whereas chi- 
ralities are of the Ising type with localized boundaries 
(cf. Fig. (2). This is also seen in the types of topologi- 
cal defects allowed by the 0(n) symmetry, namely plain 
domain walls for the Ising case, but chiral domain walls 
and vortices for XY spins^i^i^^. The classic approach of 
measuring defect energies has been to analyze the differ- 
ence between the energies Ep and Eap of systems with 
periodic (P) and antiperiodic (AP) BCs^'^. In contrast 
to a ferromagnet, however, both P and AP boundaries 
are in general frustrating for a spin glass such that, com- 
pared to the thermodynamic state, in a finite system do- 
main walls or other defects are forced into the system. 
It is therefore not quite obvious what kind of excitation 
energy the difference AE'p/ap actually corresponds to. 
For the Ising case, the difference or sum of two domain- 
wall excitations is still essentially a linear excitation of 
length L, so one might expect to observe the scaling 
|A£'|p/AP ~ L^" predicted by droplet theory^i^, which 
in fact turns out to work reasonably well for sufficiently 
large system sizes'^. The continuous case with spin and 
chiral variables is less straightforward: the mere period- 
icity of periodic or antiperiodic BCs alone might force ei- 
ther spin or chiral excitations into the system, such that 
the interpretation of the difference |Ai?|p/AP is rather 
unclear. While previously only P/AP boundaries (and 
the corresponding "reflective" boundaries for chiral de- 
fects) had been considered^i21i2^, Kosterlitz and Akinc — 
suggested to alleviate the problem of trapped defects by 
additionally optimizing over a global twist angle along 
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the boundary and comparing the resulting configuration 
to one with a relatively tt rotated or reflective boundary 
("optimum twist" method). While such extra optimiza- 
tion over a global twist minimizes the frustrating effect 
of periodicity, it docs not remove it, and hence spin and 
chiral defects remain incompletely disentangled. 

A less ambiguous setup results from comparing the un- 
restricted, open-boundary system with the case of an ex- 
plicitly inserted defect. We refer to this approach as open 
(O) and domain-wall (DW) BCa^. The DW boundaries 
could be realized by fixing the spins along one of the two 
boundary seams under consideration in the positions of 
the O ground state and the spins of the other seam in 
a configuration rotated by tt compared to the O state to 
induce a spin domain wall, cf. the sketch in Fig. [4l A 
fixing of spins is not possible within the GEM approach, 
however, since the spins do not have a direct represen- 
tation in the dual matching problem. Instead, we link 
the boundary seams periodically with strong "anisotropic 
bonds" guaranteeing a fixed relative orientation of spins. 
In the embedded matching approach, for a pair (S^, Si,) 
of spins linked across the boundary, we choose an effec- 
tive coupling 

JfiL' ~ •^strongC;^ ^1/ (8) 

instead of the general embedded coupling ([3]). Here, 
e;; = sign(5; • r), S'^ = US^. The pair (8^,8,) de- 
notes the intended final relative orientation of {Sfj,,Su), 
i.e., in terms of the angles 6, S = (cos 0, sin 0), we have 

9fj,- 9iy = 6^- dv + n for a spin DW, 

- el = 61^ + 61^ for a chiral DW. ^ ' 

Here, TZ denotes a rigid rotation of the pair (S'^,S'i/) 
such that Si, = S^. In this way, any relative orienta- 
tion {Si_i, Su) can be forced upon the system by choosing 
a (non-random) Jstrong ^ max | Jy | , such that a Jstrong 
bond is never broken in an embedded matching run (for 
our computations, we chose Jstrong = 10^). The resulting 
excess contribution of the Jstrong bonds is, of course, not 
counted in the final energy. For snapshots of excitations 
incurred by these 0/DW spin and chiral boundaries, see 
Ref. m. The prescription ® for the chiral DW bound- 
aries corresponds to a fixed axis of reflection independent 
of the disorder configuration (depending on the choice of 
origin for the angles 9i). As an alternative, we have also 
considered a minimization procedure for the choice of the 
reflection axis. In particular, we chose the reflection axis 
perpendicular to the average direction of the boundary 
spins in order to minimize the boundary perturbation. 
We found very similar results to the case of the prescrip- 
tion of Eq. ^ , such that it appears that this ambiguity 
does not have any bearing on the resulting value of the 
chiral stiffness exponent. For completeness, we also con- 
sidered the standard P/AP BCs pair as well as random 
(R) and anti-random (AR) BCs pair, where the spins are 
fixed in random relative orientations using the procedure 
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FIG. 4: (Color online) Open and domain-wall BCs. A ground 
state for open boundaries is determined, resulting in some 
configuration of the boundary spins (left). For spin domain- 
wall boundaries, the spins along the right boundary are then 
rotated by tt and all boundary spins are held fixed for a second 
ground-state computation (right). In the actual implementa- 
tion used, instead of locking the boundary spins, only their 
relative orientation is fixed. 

described above for the R computation and in relatively tt 
rotated orientations for the computations with AR BCs. 



B. Spin and chiral stiffness 

We have computed defect energies for L = 6 up to 
L = 28 systems with P/AP, R/AR and 0/DW spin and 
chiral combinations of BCs. The results are collected in 
Fig. El together with fits of the form (|A£:|)j - to the 
data. Alternativ ely, one might consid er the scaling of the 
width a{E) = ^/{{AE- {AE)j)^)j, which leads to very 
similar results, and the final estimates of this approach 
will be discussed below. For the P/AP combination, a 
pronounced crossover is observed from the sizes L < 12 
with Os = -0.724(21) to 9s = -0.433(26) for L > 16. 
L = 12 was about the maximum size considered in pre- 
vious studies^i^, and our estimate for L < 12 is indeed 
consistent with these previous results as illustrated by the 
P/AP data of Ref. 65 shown for comparison. The result 
for L > 16, on the other hand, is much smaller in mod- 
ulus and close to the value of 9s ~ —0.4 resulting from 
the "optimum twist" prescription of Ref. [H, indicating 
that it indeed reduces finite-size corrections to scaling. 
We note that the apparent crossover length L « 12 is 
comparable to the length below which no metastabil- 
ity occurs and the system behaves like a spherical spin 
glass2ii2i. The scaling of defect energies for the R/AR 
combination of boundaries behaves quite similar to the 
P / AP case (as they are still periodic in nature) , resulting 
in 9s ~ —0.519(30) for somewhat smaller systems. The 
0/DW setup, on the other hand, not hampered by the 
periodicity effect, yields an even less negative estimate 
9s — —0.207(12). Finally, for the chiral defects, we find 
9c ^ -0.090(23). 

As is seen from the variation of estimates with system 
size as well as between BCs, corrections to scaling are 
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FIG. 5: (Color online) Average defect energies {\AE\)j for 
four sets of BCs on square lattices as a function of system 
size L. The straight lines are fits of the form (|Ai5|),/ ~ 
to the data. The open black circles show the "random twist" 
result of Ref. 65 for comparison. The dashed line shows a fit 
of the form ^ with C = L> = to the P/AP data. Some 
data sets have been shifted vertically for better distinction. 



strong. A very similar situation occurs for the Ising spin 
glasa^. On general grounds^*, we expect analytic and 
non-analytic finite-size scaling corrections, 



{\AE\)j{L) = AL\l + BL-^) + C/L + D/L' 



(10) 



In contrast to the scaling ^ of the internal energy, where 
the edge and corner effects were important, the dominant 
term is here given by the non-analytic correction which 
is connected to interactions as well as self-interactions of 
domain walls^i. Due to the limited precision of the data, 
we have to set C = D = 0. We find that only for the 
P / AP case are corrections pronounced enough and a suf- 
ficient number of system sizes is available for a stable fit 
that includes a non-analytic correction term. For 
this case (including all data points L ~ 6, . . . , 28), we ar- 
rive at 0s = -0.18(27), LJ ^ 1.21(28) with Q = 0.43 (cf. 
the dashed black curve in Fig. [5]) , whereas a fit includ- 
ing all data without the correction term (i.e., addition- 
ally setting B — 0) yields an unacceptable quality-of-fit 
Q = 6 X 10^^^. Hence, including corrections to scal- 
ing seems, in principle, able to reconcile the data from 
different sets of boundaries. However, this can only be 
done at the expense of mounting statistical uncertainties 
which cannot be afforded for useful estimates with the 
available level of accuracy and range of system sizes. 

As an alternative systematic treatment of scaling cor- 
rections for domain-wall calculations, it has been sug- 
gested to consider L x M systems with variable aspect 
ratios R = M / L (with the change of BCs occurring along 
the edges of length L)^. In general, scaling depends on 
the aspect ratic''-, so we can write 



\AE\)j{L)=L<'F{R), 



(11) 



with some scaling function F{R). In the extreme case 
of very elongated systems, R — > oo, one should approach 
the behavior of the one-dimensional (ID) system. For 
the ID Ising model, it is found that (|A£;|),7(L) - 1/M 



independent of BCgdS. For the XY symmetry, the situa- 
tion is more involved due to the simultaneous presence of 
spin and chiral variables. It appears that there are dis- 
tinct exponents 9 — (d — 2) = — 1 related to spin waves 
and 6 w —1.9 related to chiral excitations for the (quasi) 
one-dimensional XY spin glass^SiS^. We are not aware, 
however, of any tests of universality with respect to dif- 
ferent BCs probing the same excitations in this model. 
However, it is clear that the influence of the boundaries 
must diminish for i? ^ 1 as the average distance to the 
boundaries with variable conditions increases. To test 
this hypothesis, we performed additional GEM compu- 
tations for systems with aspect ratios R = 2 and R = 6 
with sizes L = 4, 6, . . . , 16 (i? = 2) resp. L = 3, 4, . . . , 9 
{R = 6) and 5 000 disorder replica per lattice. In Fig. [H] 
we present the results for domain- wall energies as a func- 
tion of system size L at aspect ratio R — 6. As is appar- 
ent, the slopes of the data for boundaries probing spin 
excitations, i.e., P/AP, R/AR and 0/DW spin boundary 
conditions, are quite compatible with each other asymp- 
totically, whereas energies of chiral domain walls decay 
considerably slower. If, as it hence appears, there is in- 
dependence of BCs as i? — > oo, one expects scaling cor- 
rections depending on boundaries to gradually disappear 
as more and more elongated systems are being consid- 
ered. For the Ising case, it was observed that to good 
approximation this limit is attained as^i^^ 



0{R) = 0{R = oo) + Ar/R. 



(12) 



Our estimates of 9{R) for the data at aspect ratios R = 1, 
2 and 6 for the different sets of boundary conditions are 
presented in Fig. [7] together with fits of the form (fT2)) 
to the data, also including the alternative scaling of the 
width (j{E). For each aspect ratio, corrections at fixed R 
were taken into account by omission of points from the 
small-L side for the fits of the form (fTU|) instead of includ- 
ing any correction terms (i.e., we set _B = C = 13 = 0), 
cf. the fits presented in Fig. [Sj In particular, data points 
were successively omitted from the small-L side until a 
satisfactory quality-of-fit Q was achieved. We indeed find 
a clear convergence of the P/AP, R/AR and O/DW spin 
stiffness exponents as more elongated systems are con- 
sidered, indicating independence of BCs in this limit. 
We collect the extrapolated estimates in Table HI finding 
them all well compatible within statistical errors. As our 
final estimate we cite the weighted mean 0s = —0.329(14) 
(we do not include the a{E) results in the average since 
they are not statistically independent). The chiral stiff- 
ness exponent 6c, on the other hand, seems to be al- 
most independent of the aspect ratio (at least as far as 
the scaling of the difference |AS| is concerned) with an 
asymptotic value of 0c{R = oo) = —0.114(16), clearly 
different from the value found for the spin stiffness, indi- 
cating spin-chirality decoupling in this model. 
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\AE\ 


a{E) 


P/AP 


-0.338(20) 


-0.348(19) 


R/AR 


-0.329(28) 


-0.346(13) 


O/DW spin 


-0.308(30) 


-0.306(26) 


O/DW chiral 


-0.114(16) 


-0.095(38) 



TABLE I: Spin and chiral stiffness exponents d{R = oo) from 
different sets of BCs as extrapolated via aspect-ratio scaling 
according to (fT2)) for |A_E| and the width u{E). 



FIG. 6: (Color online) Scaling of average defect energies 
(|A_E|)j for four sets of BCs at aspect ratio i? = 6. The 
straight lines are fits of the form (|A_E|)j ^ to the data. 
Data sets have been shifted vertically for better distinction. 
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FIG. 7: (Color online) Dependence of the stiffness exponents 
9s and 9c on the aspect ratio R of the systems for different 
pairs of BCs. The solid lines show fits of the form (|12|l . 



C. Domain-wall fractal dimension 

On general grounds, the domain walls or droplet sur- 
faces of the Ising spin glass have been assumed to be me- 
andering curves with a non-trivial fractal dimension ds 
(see Refs. [Uli). Measurements for the case of the 2D 
Ising spin glass with Gaussian interactions indeed yield 
a non-generic value^l ds « 1.27, and for the bimodal dis- 
tribution quite similar estimates are foun d^^i^^i^°° . A 
similar value has also been found in two-dimensional 
0{n) spin glasses with additional random (off-diagonal) 
anisotropyiSi, which are expected to belong to the Ising 
spin-glass universality classiS^. The fractal dimension of 
domain walls is relevant to the question of whether the 
mean-field scenario, which predicts space-filling domain 
walls, i.e., dg = d, applies to low dimensions'^. Consider 
the link overlap. 



1 

iVi 



{Sf-Sf)iSf 



sf), 



(13) 



where A^i denotes the number of bonds of the lattice, 
and the superscripts A and B stand for two replica of 
the system, here for the two BCs for a system with a 



given realization of the random bonds. For the Ising 
case, bonds crossed by the domain wall contribute —1 to 
qf^ and all the others -1-1, such that 



(1 



1l I J 



(14) 



which allow s to determine dg numericalljs^l. The authors 
of Ref. Il03l used the same relation for continuous spins 
(extrapolating from finite temperatures to T = 0). Due 
to the continuous nature of the spin response to the per- 
turbation on the boundary, however, all spins contribute 
here to the reduction of qf^ and it is therefore almost 
inevitable that (1 — qf^).j ^ L^, implying ds — d. (It 
is likely that the small difference d — dg found by the 
authors of Ref. ilOS is an artefact of the extrapolation 
from finite temperatures.) Here, we indeed find, e.g., 
d - ds ^ 0.0010(14) for the P/AP case and similar re- 
sults for the other boundary conditions. Due to the above 
argument, however, this should not be seen as evidence 
for space-filling domain walls, but rather as an indicator 
for the inadequacy of the method. 

As discussed above, domain walls in the system con- 
sist of the union of boundaries of rigidly 0(n) rotated 
clusters. To reveal this structure in the ground states 
identified, firstly one should get rid of the global symme- 
try by rigidly 0(n) rotating one configuration such as to 
maximize the total overlap q^^ ~ ^ Qap^afi^ where 
denotes the corresponding global rotation matrix, 

and 



1 



(15) 



is the matrix of (site) overlaps. This is done by a singular 
value decomposition to diagonalize q^p , in which case 
q^^ is just the trace of the resulting diagonal q^^ , see 
Ref. m. We refer to q"^^ as the "optimized" overlap (or 
projection) in contrast to the plain scalar overlap q^^ — 
X^a laa ■ Ouce this global symmetry has been removed, 
one can attempt to determine local rotation matrices to 
reveal rigid relative rotations of clusters of spins between 
the pair of configurations considered. To this end, we 
consider a locally averaged overlap matrix, 

q^^ix) = Z{x)-^ "^(^^ - ^)StaS?P^ (16) 

i 

where Z{x) = wixi — a;), and w{xi — x) is a rapidly 
decaying weighting function. The optimal local rotation 
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FIG. 8: (Color online) Scaling of the chiral link overlap qi ^ 
of Eq. Hll and the spin overlap qfj^ of Eq. HH]) for L x L 
P/AP systems. The lines show fits of the functional form (I14p 
to the data. 

matrix R^^ (x) again follows from a singular value de- 
composition. The averaging is necessary here since the 
comparison of a single pair of spins does not allow to 
uniquely determine an 0(2) matrix. In particular, one 
could not decide whether a proper or an improper ro- 
tation is involved. The weighting function wlxi — x) is 
chosen exponentially decaying and cut off after two or 
three lattice spacings. The following results are indepen- 
dent of any fine-tuning of this falloff. The snapshots of 
Fig. [2] show the angles and signs of the thus determined 
rotation matrices R^p{x) (for snapshots of domain walls 
induced by a change of BCs, see also Ref. losh . Domain 
boundaries can be defined from R^^{x) separately for 
the spin and chiral parts: for the latter any bond with 
different signs of the determinant of R'^p{x) at its two 
ends is crossed by a chiral domain wall, and for the for- 
mer any bond with a change of the rotation angle (p of 
R^p{x) by more than a threshold value is crossed by 
a spin domain wall. To determine for the chiral case, 
we define a link overlap analogous to p3|) . 

itx = ]^ J2^detR^^{x,)] [detR^^{x,)l (17) 

measuring the change in rotation signs along bonds. An 
analogous expression for qf^ is used for df of the spin 
domain walls, 

C = E^*o(^.,^.), (18) 

where 9^g{xi,Xj) = — 1 if the angle between Si and 
Sj exceeds 0o and detil'^^(a;i) = det i?^'^(a;j), and 
0(j,^{xi,Xj) = otherwise. Hence, bonds including a 
sign change of AetR^^{x) are explicitly excluded from 
the counting to optimally decouple spin and chiral mea- 
surements. Fits of the form ([H]) are used to determine 
the corresponding fractal dimensions. 

We consider the same sets of BCs and system sizes used 
for the determination of the stiffness exponents above. 



Figure [8] shows the L dependence of the spin and chi- 
ral link overlaps qfj^ and qf^ ^ respectively, for the case 
of P/AP BCs and L x L systems, revealing clear scal- 
ing signals. For the chiral overlap, a fit of the form (fH|) 
to the data yields = 1.4291(44). The 0/DW chi- 
ral boundaries presumably should show the clearest sig- 
nal since direct visual inspection of the configurations 
reveals that typically a single domain wall roughly paral- 
lel to the variable boundary is excited. For this case we 
arrive at d^ = 1.394(12). The sizable, although not dra- 
matic, difference for the two choices of BCs shows that, 
as expected, still only partially controlled finite-size ef- 
fects are present. Unfortunately, aspect-ratio scaling is 
not a suitable technique for the determination of d^, at 
least for the system sizes accessible here: for the sizes 
L = 3, . . . , 9 of the R — 6 systems, the wall simply has 
no space to properly "meander" over the length of this 
couple of lattice spacings, resulting in very strong finite- 
size effects in L. Finally, as the third combination of 
boundaries considered, the 0/DW spin boundary setup 
yields an extremely weak signal for the chiral domain 
walls since improperly 0{n) rotated domains occur only 
scarcely there, and we hence cannot produce a reliable 
estimate. We take the average of the P/AP and 0/DW 
chiral results with sufficiently strong signals as our fi- 
nal estimate, = 1.425(12). The analysis of spin do- 
main walls from the overlap qf^ analogous to Eq. (fT7)) 
is more involved due to the general dependence of the 
results on the cutoff angle if it is chosen too small, 
smooth rotations are counted as parts of domain bound- 
aries, ultimately leading to space-filling domain walls in 
the limit (/)o 0. Choosing it too large, on the other 
hand, misses to recognize some of the cluster bound- 
aries. We find somewhat stable results, however, in a 
range 0.5 < (j)o < 0.75, and estimate = 1.2379(97) for 
P/AP and df = 1.251(21) for 0/DW spin boundaries. 
The 0/DW chiral conditions yield too weak a signal for 
the determination of df, and we again cite the average 
(if — 1.240(21) of the remaining two sets of boundaries 
as our best estimate. Thus, to the extent that finite-size 
corrections and the cutoff dependence of df are under 
control, it appears that spin and chiral domain walls have 
different fractal properties. 



V. OVERLAPS 

One of the most striking results of the mean-field the- 
ory of the SK model is the occurrence of a multitude of 
pure states, resulting in a non-trivial distribution of the 
(thermally and disorder) averaged overlap parameter—. 
Consequently, estimates of the form of the overlap distri- 
bution have been regularly considered as crucial bench- 
marks of the extent to which replica symmetry break- 
ing occurs in finite-dimensional spin glasses^. In recent 
years it has become clear, however, that great care has to 
be taken when drawing conclusions from measurements 
on finite systems for the structure of pure states in spin 
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glasse a^^i^^'^"^'^"^ . According to basic principles of sta- 
tistical mechanicsiS^, a thermodynamic state should be 
defined as a limiting probability measure from a series 
of systems of increasing size, where a specific state is 
selected by an accompanying series of BCs. The result- 
ing state is then identified by the values in this limit (if 
they exist) of all possible correlation functions in a cen- 
tral volume of the system, far away from the perturbing 
presence of the boundaries. If such a state cannot be 
decomposed into other thermodynamic states, i.e., it is 
mixed, it is termed pure. While this is rather clear in ho- 
mogeneous systems such as ferromagnets, unfortunately, 
it is not known how a pure state should be explicitly con- 
structed in spin-glass systems^i. In particular, there is 
"chaotic size dependence" in such a sequence of system 
sizes resulting in oscillating and non-convergent correla- 
tion functions. It is also unknown how a convergent sub- 
sequence could be selected by coupling-independent BCs. 
It is clear, then, that considering the overlap distribution 
in whole, finite systems is not an appropriate way to de- 
tect the presence of a multitude of pure statesiiS. In par- 
ticular, the presence of domain walls separating patches 
of pure state configurations in finite systems might lead 
to a non-trivial overlap distribution even if only a single 
class of pure states related by a global 0(n) transfor- 
mation was to exist. If, on the other hand, spins in a 
volume far away from the boundaries always converge 
to the same configuration (up to a common rotation) as 
L — > oo, independent of BCs, a multitude of nonequiva- 
lent and observable pure states would appear to be ruled 
out. 

An asymptotic independence from BCs of spin con- 
figurations in central windows can be detected either 
in (window) overlap distributions or in the behavior of 
correlation functions. These approaches are essentially 
equivalenli^. To see this, consider two spins 5"^ and Sj 
inside a window W containing Nw = spins (for sim- 
plicity we consider Ising spins first) . The average squared 
difference in the two-point functions between two systems 
with different BCs labelled A and B, 

^ = ^T. - {S^Sf)Tf)j, (19) 

by construction scales to zero at fixed window size w if 
and only if the two-point functions for all window spins 
agree for almost-every disorder configuration (here {■)t 
denotes a thermal and (•) j a disorder average). It is then 
easy to see that in terms of the optimized overlap q"^^ 
of window spins (which for the Ising case n = 1 is just 
|(7"^"^|) this is equivalently written as 

^^m^^f + {q^''f-2{q^^f)T,j. (20) 

For the case of our XY model with a unique state at 
T = 0, up to a global 0(n) rotation (see Sec. lIII Bl above). 
q"^^ and q^^ are of course always unity. Higher order 
(even) correlation functions lead to similar expressions 
with higher moments of q"^^ , such that in the most gen- 
eral context, the vanishing of the difference /S.P^^{q) = 
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FIG. 9: (Color online) Sampled distribution of optimized 
overlaps q"^^ in central windows of size w = 4 for a change 
of BCs from periodic to antiperiodic. As the system size is 
increased, the distribution accumulates all its weight close to 



P^^ (q) + P^^ (q) - 2P^^ (q) of overlap distribution fimc- 
tions is equivalent to the invariance of window correla- 
tion functions under a change of BCs. For the present 
case of non-degenerate ground states, this vanishing is 
attained only for the trivial form P^^{q) = 5{q — 1), i.e., 
if and only if the window configurations are completely 
unchanged for almost every disorder configuration. A 
calculation along the same lines establishes the analo- 
gous relation between correlation functions and overlaps 
for the case of continuous spins. 

We did not attempt to sample "all" BCs, but used the 
P/AP, R/AR and 0/DW boundaries already considered 
for the defect-energy calculations. If the size w of the 
overlap window comes too close to the system size L, 
a crossover to the behavior of whole-system overlaps is 
seen, such that one strives to guarantee w <C i. For too 
small windows, on the other hand, discretization or finite- 
size effects appear, and there is a certain probability that 
the two window configurations can be made to (almost) 
perfectly overlap by an 0(n) rotation although they are 
affected by the influence of the changed boundaries — 
down to the extreme case of w = 1 where each pair of 
configurations has q^^ = 1. Thus, the intended behavior 
should appear for window sizes 1 <^ w <^ L. As we shall 
show below, with available system sizes, these restrictions 
are best fulfilled for about 4 < w < 8, but finite-size ef- 
fects are still clearly visible due to the maximum L « 28 
that we can reach with the GEM. In Fig. [9] we present the 
sampled distribution function of window overlaps q'^^ 
for the P/AP combination of BCs with w = 4, showing 
a clear tendency of statistical weight to accumulate at 
q^^ = 1 as the system size L is increased, suggesting a 
convergence to the trivial form P^^{q) = 5{q — 1) in the 
limit L —f oo. The overlap distributions for the remain- 
ing pairs of BCs show a qualitatively similar behavior. 
For a more systematic study of this apparent conver- 
gence, we consider the total probability p'^'^(L, w) for the 
occurrence of an overlap q^^ smaller than a threshold 
value qo, measuring the distribution weight away from 
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FIG. 10: (Color online) Scaling of the probability p^^{L,w) 
of an overlap q^'^ < 0.9 for the P/AP pair of BCs as a function 
of L/w. The lines show fits of the form pil) at fixed w to the 
data. 

q-^^ = 1. To the extent to which the configurational 
changes induced by a change of BCs take the form of 
domain walls, this just measures the probability of the 
window being crossed by one of them. As discussed in 
Sec. IIV C\ such domain walls are fractal curves occupy- 
ing a volume ^ such that a fraction of the lattice 
proportional to will be touched by a domain 

wall (whether the relevant fractal dimension would rather 
be connected with spin or chiral defects will have to be 
decided post factum, see the discussion below). Since 
the relevant length scale in fact is L/w (remember that 
fractals are scale invariant), one expects an asymptotic 
scalingiSi 

p^^(L,H=Qo(LM"^'"'=^- (21) 

Figure Uni shows the scaling of the crossing probability 
p^^{L,w) for go = 0.9 and the P/AP pair of BCs. We 
additionally considered cutoffs go — 0.95 and go — 0.99, 
changing the prefactor Qo of ([2T|) and slightly modifying 
corrections, but yielding the same leading (L/k;)"'^''^''^^ 
behavior. There is a clear signal of scaling of p^^{L, w) 
to zero as the system size L is increased at fixed window 
size w. The expected collapse of data for different win- 
dow sizes, however, only appears at u; > 8, at the edge 
of the window sizes that can be reasonably considered 
with the system sizes at hand. From fits to the scaling 
form ([^T|) we find exponent estimates ds — 1.442(29), 
ds = 1.444(16) and 4 = 1.472(21) for window sizes 
w — A, w = 6 and w — 8, respectively. These are 
not only consistent with each other, but also compati- 
ble with the estimate d} — 1.425(12) found above for the 
fractal dimension of chiral domain walls. The remain- 
ing combinations of BCs yield rather similar results, with 
somewhat more pronounced finite-size corrections for the 
L < 20 systems considered there (in contrast to the 
aspect-ratio scaling of the domain-wall energies, where 
P/AP boundaries entailed the stronger corrections). In 
particular, we find dg = 1.404(17) from R/AR boundaries 
(w = 6), 4 = 1.491(28) from 0/DW spin BCs {w = 4) 
and ds = 1.489(59) for the 0/DW chiral combination 




FIG. 11: (Color online) Spatial distribution of the internal 
field per link in the ground state of an example 28 x 28 disorder 
configuration with open-periodic BCs. There is a clear inverse 
correlation with the density of frustrated plaquettes (small 
white squares). 

{w — 6). It thus appears that the relevant domain walls 
seen in the window overlaps are of the chiral type. In 
total, clear signals of a scaling of the weight away from 
q^^ = 1 to zero are seen in all cases, indicating a "de- 
flection to infinity"— of domain walls in the center of 
sufficiently large samples. The estimates for the fractal 
dimension seem to preclude an excitation of space-filling 
domain walls (i.e., ds = 2) as predicted by a many-state 
picture at least from the coupling-independent and sym- 
metric BCs considered. This is in agreemen t wit h the 
findings for the 2D Ising spin glass of Refs. [5llll04 

VI. MAGNETIC PROPERTIES 

A. Internal Fields 

It is straightforwardly seen that a necessary condition 
for the Hamiltonian ^ to be in a metastable state is 
that each spin is parallel to its local internal field 

f< = -J2j^i^i' (22) 

which for convenience has been normalized to represent 
the field per link by dividing by the local co-ordination 
number qi, accounting for the difference in connectivity 
between spins along open boundaries and in the bulk. 
An iterative alignment of spin orientations along the in- 
ternal field directions indeed drives the system towards 
a metastable (but not a ground) state, and the corre- 
sponding minimization technique is referred to as the 
Walker- Walstedt or spin-quench algorithmlS. The local 
strength of internal fields is closely linked to the presence 
of frustration in the system, with a weak internal field 
indicating a local uncertainty of the system with respect 
to choosing a given spin's orientation such as to mini- 
mize the total energy. We indeed find a strong inverse 
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correlation between the size \h\\ of the internal fields in 
the ground state and the density of frustrated plaquettes 
of the system. This is illustrated for the ground state 
of an example disorder configuration of 28 x 28 spins in 
Fig. 111! It is clear that in terms of finding a ground 
state, optimal spin orientations in regions of low frustra- 
tion (and consequently large internal fields) are rather 
easily determined, while it is much harder to optimize 
highly frustrated areas. This connection was exploited in 
the improved spin-quench method suggested in Ref. [3^ . 
One then also expects a certain correlation of the bound- 
aries of rigid clusters discussed above in Sec. Illll and the 
occurrence of a high density of frustrated plaquettes. 

Speculations on the occurrence of a certain propor- 
tion of sites with almost or exactly vanishing internal 
fields played some role in early models of the specific 
heat in spin glasses, see Ref. UQ and references therein. 
For the infinite-range case, it was subsequently estab- 
lished that \h\\ > hg > 0, such that (almost) free spins 
do not occur-i2i. For the short-range Heise nberg rn odel, 
a similar exact bound was reported in Refs. Il08lll09l . For 
the bimodal Ising spin glass, on the other hand, a mul- 
titude of flippable spins exists in each ground stateii^. 
For the intermediate XY case, however, to our knowl- 
edge no exact results are available. Figure [12] shows the 
distribution of internal fields for 28 x 28 systems. In- 
dependent of BCs, no spins with fields < 0.6 = 
occur. Smaller systems have slightly larger gaps h^, but 
the size dependence is very weak such that it is clear that 
a finite /ig survives in the thermodynamic limit. The ab- 
sence of almost free spins in the ground state clearly il- 
lustrates the fact that there are no localized low-energy 
excitations and, instead, the lowest excitations are the 
extended defects corresponding to domain walls or 0(n) 
rotations of rigid clusters, as discussed in Sections llIII and 
IIVI The distributions P(|/j.'|) for the ground state show 
a characteristic behavior with the bulk of spins having 
1.25 < \h\\ < 1.75 indicating strong frustration, but a 
certain number of almost fully satisfied spins leading to 
an extra peak close to \h\ \ « 2 associated to regions free 
of frustrated plaquettes. There are only small differences 
between the systems with open and periodic BCs, most 
noticeably a relative enlargement of this peak of unfrus- 
trated spins for the open boundaries resulting from the 
relatively lower frustration of the extra boundary spins. 
For comparison, we also show the corresponding distri- 
bution for states found from the Walker- Walstedt zero- 
temperature quench, which asymptotically should result 
in an even sampling of all metastable states. This leads 
to a rather different distribution of internal fields where 
much less of the inherent frustration has been relieved. 
These curves for the finite-range model are markedly dif- 
ferent from the distribution (again over all metastable 
states) of the infinite-range 0(2) SK modeli^l shown for 
comparison which, due to the lack of local fiuctuations in 
this infinite-dimensional model, is much more symmetric 
and more evenly frustrated. 
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FIG. 12: (Color online) Average distribution of internal 
fields The "GEM" curves represent ground states of 

28 X 28 systems with open-open and open-periodic BCs. 
The "quench" data correspond to the average over a broad 
range of metastable states as found from the Walker- Walstedt 
algorithm/'^. The curve labeled "SK model" shows the corre- 
sponding distribution in the 0(2) SK model^. 



B. Magnetic phase diagram 

For small concentrations x and/or rather weak im- 
purity bonds, one might expect the system to exhibit 
some rather extended short-range order—, or perhaps 
even long-range order~, and from example calculations 
we indeed find clearly non-zero magnetizations for this 
region, whereas strongly disordered states of spin-glass 
signature appear once a larger number of strong impu- 
rity bonds is being introduced. Analytical calculations 
regarding the question of order and canting induced by 
impurity bonds have only been possible for isolated de- 
fects or using the mean-field type coherent potential ap- 
proximation (CPA) of unclear precision for the case at 
hand^Si^i^. As a preliminary exploration of the zero- 
temperature phase diagram, we consider the modulus of 
the magnetization, {\m\)j, where m — Si)/NQ. Fig- 
ure [T3] shows a contour plot of (|m|)j as a function of 
X and A [cf. Eq. (|4])] as found from GEM ground-state 
computations for systems of 16 x 16 spins and open BCs. 
Numerical calculations were actually performed at the 
grid of 30 points indicated by the tick marks in the plot, 
and the contours were produced from a S"^*^ order poly- 
nomial interpolation of the results. The bold black line 
indicates the locus of the curve {\m\)j — 0.5 which could 
be used as a possible definition of a finite-size approxi- 
mation of the phase transition line, approaching the true 
transition line as L ^ oo. The prediction of the CPA, 
namely^ 



1(WA)^ 
"^^^^^-2 1 + A ■ 



(23) 



is indicated by the bold white line. In the limit A —> 
of a diluted and unfrustrated system, order must disap- 
pear at the (bond) percolation threshold Xc{X = 0) = 0.5 
which is (by chance) already very well reproduced by the 
chosen definition (|m|)j = 0.5 at the system size dis- 
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FIG. 13: Contour plot of the average ground-state magneti- 
zation (|m|)j for 16 X 16 systems with open BCs as a function 
of a; and A. Darker shades indicate larger magnetizations, and 
contours are separated by Sm = 0.1. The bold black line de- 
notes the locus of the curve (|m[) j = 0.5, and the bold white 
line corresponds to the CPA prediction of Eq. (|23|l . 



played. Also, Xc(0) = 1/2 for the CPA form There 
is no agreement in the literature about the general form 
of Xc{X) as a: — *■ 0, i.e., in the experimentally most rel- 
evant regime of small dilution. According to the 
CPA predicts a finite intercept Ac = 1 for x = 0, and this 
claim was supported by the interpretation of the average 
magnetization found from a local spin-quench computa- 
tion in Ref. m (at the same time, however, further results 
from the same study seemed to indicate the absence of 
long-range order for the whole phase diagram). An RG 
calculation for the non-linear a model'^'^, on the other 
hand, appeared to indicate Ac oo as a; — )■ 0. The 
results of Fig. [13] certainly rather suggest that Ac > 1 
for a; ^ 0, but this impression might be distorted by 
finite-size effects. The question of long-range order is 
here rather subtle"^^, and it certainly cannot be excluded 
from the data shown that the indicated transition line is 
related to a crossover rather than a phase transition. In 
any case, this picture, of course, needs to change again 
at finite temperatures, where true long-range order is im- 
possible due to the continuous symmetry of the system. 
A more thorough, finite-size scaling analysis of this ques- 
tion, together with an investigation into the nature of the 
apparently ordered phase (i.e., whether it is accompanied 
by transverse spin-glass order'^^), is left as an interesting 
future application of the GEM technique. 



VII. CONCLUSIONS 

We have shown that a newly developed specially tai- 
lored heuristic optimization technique combining exact 
ground state computations for embedded Ising variables 
with a genetic algorithm allows for a thorough investiga- 
tion of the ground state of the planar EA spin glass in 
two dimensions, contributing to a resolution of a number 
of long-standing issues for this system. 



The main result consists of a careful determination 
of the spin and chiral stiffness exponents of the model 
with bimodal coupling distribution, which characterize 
its critical behavior and directly relate to the question 
of a possible spin-chirality decoupling in this system. In 
principle, a stiffness exponent related to domain walls 
can be extracted from scaling corrections to the finite- 
size behavior of the ground-state energy, but the preci- 
sion of the estimate resulting from this method is not 
satisfactory. Instead, we investigated the properties of 
domain-wall defects directly induced by a suitable change 
of boundary conditions. In response to doubts regarding 
the suitability of the conventionally considered combi- 
nation of periodic and antiperiodic boundary conditions 
resulting in an effective measurement of differences of 
defect energies^^, we explicitly excite domain- wall de- 
fects of the spin and chiral type relative to a system with 
free boundaries by applying special "domain- wall bound- 
ary conditions" . Corrections to scaling for these as well 
as for the additionally considered periodic-antiperiodic 
and random- antirandom boundaries are found to be large 
and their resolution, from explicitly including correction 
terms in fits to the defect energies, is not fully satisfactory 
given the accessible system sizes. However, application 
of the aspect-ratio scaling technique^^ allows for a good 
control over the boundary-dependent corrections, lead- 
ing to a consistent estimate Og — —0.329(14) of the spin 
stiffness exponent for the different sets of boundary con- 
ditions used. This value is considerably smaller in mod- 
ulus than the results of previous authors, quoting values 
of Bs = -l.O.S"* -0.9^ -0.8^ -0.78^ and -0.37(2)^ 
respectively. The spread of these previous results nicely 
illustrates the magnitude of present finite-size corrections 
as well as the hardness of computing true ground states 
for sufficiently large systems. The fact that in the present 
work, consistency between estimates of the stiffness expo- 
nent from different sets of boundary conditions with cor- 
rection amplitudes of opposite signs could be achieved, 
makes it appear plausible that finally scaling corrections 
are under control to an acceptable degree of p recision. 
Concerning the claim put forward in Ref. |65| that the 
P/AP combination of boundary conditions cannot be 
used to estimate the spin stiffness exponent, the results of 
the aspect-ratio scaling show that the assumed trapping 
of domain walls in periodic systems^^ only affects scaling 
corrections, but does not alter the asymptotic stiffness. 
This value of Og is remarkably close, but apparently not 
identical to the estimate 6s — —0.287(4) found for the 
two-dimensional Ising spin glass with Gaussian coupling 
distribution^''. Also, scaling corrections to the ground 
state energy are found to be very close to the Gaussian 
Ising model, but rather different from the bimodal Ising 
spin glass. Chiral boundaries, on the other hand, lead 
to a clearly different exponent Oc = —0.114(16), which 
is again remarkably smaller in modulus than previous 
estimates of Oc = -0.6fiii -Q.ZS^^ and -0.37(1),^ 
respectively. Consequently, it appears that spin and chi- 
ral degrees-of-freedom order with different exponents in 
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this system although, of course, a study for finite system 
sizes can never guarantee to probe the true asymptotic 
behavior. Naively, one might expect that the Ising-like 
chiral variables should lead to the value 6'c = of the 
bimodal Ising spin glass^^. It should be noted, however, 
that the chiral variables interact with an effective long- 
range Coulombic interaction^^, possibly leading to a shift 
in universality class. On the other hand, for the recently 
established spin-chirality decoupling in the periodically 
frustrated XY model^i^, the chiral transition appears 
to show Onsager exponents, but this behavior is only 
seen for huge lattice sizes of L « 10"^, with wildly differ- 
ent behavior in the crossover region^. For our data, since 
for a 9c very close to zero it becomes very hard numeri- 
cally to distinguish the leading scaling from sub-leading 
corrections, we cannot definitely exclude the possibility 
of 9c — 0. In fact, including a finite asymptotic value 
AEq in the form pUj) fits consistent with 9c — can be 
achieved as well. In the context of the observed spin- 
chirality decoupling, it is therefore quite interesting that 
we find a 9s value so close to that of the 2D Ising spin 
glass with Gaussian interactions. 

We find that the structure of metastable states of the 
model can be understood in terms of metastates of clus- 
ters of rigidly locked spins, which undergo common rela- 
tive 0(n) rotations between different metastable states. 
To each metastable state belongs a spectrum of contin- 
uous spin- wave excitations leaving the associated energy 
minimum. These will be considered in a future inves- 
tigation. Low-energy excitations within the manifold of 
metastable states, on the other hand, are extended ob- 
jects, corresponding to 0(n) rotations of one or several 
rigid clusters. Almost free spins corresponding to local 
excitations do not occur. Many such metastable configu- 
rations are found closely above the ground state, leading 
to a quasi-continuous spectrum in the thermodynamic 
limit. The ground states themselves, however, are found 
to be unique up to a global 0(n) transformation, imply- 
ing a critical exponent = of the ground-state corre- 
lation function. Together with the estimates for 9s and 
9c, this gives a rather complete characterization of the 
T = critical point of the model. 

We investigated the nature of elementary excitations 
induced by a change of boundary conditions by analyzing 
the locally averaged optimal rotation matrices between 
ground-state configurations. Chiral domain walls are by 
definition sharply localized. We find them to be fractal 
curves, and from the scaling of the corresponding link 
overlaps we estimate a fractal dimension d} — 1.425(12). 
The definition of spin domain walls depends on a some- 
what arbitrary cutoff angle, but we find the fractal di- 
mension df = 1.240(21) to be independent of the cutoff 
within some reasonable range. Recently, evidence was 
presented for the hypothesis that domain walls in the 
Ising spin glass are special fractal curves described by 
"stochastic Loewner evolution" (SLE)ii2iiii^. Compar- 



ing the numerical results with conformal weights of the 
Kac table, the authors of Ref. Ill2l conjecture a relation 
ds = 1 + 3/4(3 + 9). This we find fulfilled within sta- 
tistical errors for df, but not for the chiral exponent 
where one might have rather expected it to hold. In any 
case, from the excitation via a change of boundary condi- 
tions, we do not find space-filling domain walls {dg — 2) 
one necessarily expects for a model with many thermo- 
dynamic pure states2i. This appears to be corroborated 
by an explicit investigation of the distribution of over- 
laps in central windows far away from the boundaries, 
which, for the system sizes accessible here, appears to 
converge to a trivial form P^^{q) — 5[q — 1). This 
shows that the states defined by the considered set of 
boundary conditions always converge to the same class 
of thermodynamic pure states related by a global 0(n) 
rotation, and we do not find evidence for a multitude of 
pure states. Naturally, however, consideration of differ- 
ent types of boundary conditions might yield different 
results. The rate of this convergence is again related to a 
domain-wall fractal dimension, with a value compatible 
with c?J as determined directly. 

A multitude of extensions of the present work comes 
to mind. Besides the mentioned investigation of spin- 
wave excitations complementing the information on 
metastable states presented here, certainly the magnetic 
phase diagram discussed in Sec. IVI[ which is relevant to 
a number of experimentally realized systems, deserves 
a more extensive investigation. Concerning the struc- 
ture of the ground and excited states in terms of clusters 
of rigid spins, a more systematic anal ysis of the mor- 
phology along the lines of Refs. [85llll4l appears promis- 
ing. Clearly, investigations of Gaussian bond distribu- 
tions and of the Heisenberg case would be highly interest- 
ing enterprises. Above all, however, the important ques- 
tion of spin-chirality decoupling in three-dimensional spin 
glasses could be tackled with a suitable generalization of 
the GEM approach employed here, replacing the embed- 
ded matching component inside the genetic algorithm by 
a technique applicable to three-dimensional systems. 
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